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ABSTRACT 



Time evolution of quantum tunneling is studied when the tunneling system is 
immersed in thermal medium. We analyze in detail the behavior of the system after 
integrating out the environment. Exact result for the inverted harmonic oscillator 
of the tunneling potential is derived and the barrier penetration factor is explicitly 
worked out as a function of time. Quantum mechanical formula without environment 
is modifed both by the potential renormalization effect and by a dynamical factor 
which may appreciably differ from the previously obtained one in the time range of 
l/(curvature at the top of potential barrier). 
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I Introduction 



Tunneling phenomena in thermal medium are of both theoretical and practical 
interest in many areas of physics. For instance, in cosmology there are a variety 
of tunneling phenomena that may have occured in the evolution of our universe. 
If the phase transition of the electroweak gauge symmetry is of the first order as 
assumed in the electroweak scenario of baryogenesis Jl|, then the tunneling from a 
metastable state to the true ground state of the electroweak theory must go through 
either quantum effect or thermally activated barrier crossing. Another possibility of 
the first order phase transition in cosmology is the quark-hadron phase transition. 
We should neither fail to mention the classic example of tunneling phenomenon that 
takes place in the central core of stars: thermonuclear reaction M. 

A common problem to all these is the presence of environment: the tunneling 
we are interested in does not take place for the system in isolation. Under this cir- 
cumstance we are very much interested in how substantially quantum mechanical 
formula for the tunneling rate is modified by dissipative interaction with surround- 
ing medium. Despite of this obvious interest, many past works in cosmology and 
astropysics have relied on simple methods to deal with the tunneling problem, either 
by using the bounce solution in the Euclidean approach [[|, || or by exploiting 
some variant of the classical Vant Hoff-Arrhenius law 0,0. 

On the other hand, since the pioneering work || on quantum dissipation to tun- 
neling phenomena at zero temperature, there have appeared many extended works 
in condensed matter or statistical physics community (a partial list of these works 
is given by || , [Kj , [|ll|] ) . Intensive theoretical activities in this field are presumably 



related to experimental possibility of observing the macroscopic quantum tunneling 
in various areas of condensed matter physics. The problem is however not simple 
and only a limited class of problems have been addressed. Thus even in idealized 
models one often assumes that the entire system is in thermal equilibrium and at- 
tempts to derive quantities of limited value such as the decay rate of the metastable 
state by an extention of the bounce analysis. In some of these works M, |TU|, the key 
quantity is the imaginary part of the free energy, which may be interpreted as the 
decay rate. However, since equilibrium value of the free energy of the entire system 
is necessarily real, one must extract the imaginary part by projecting to the initial 
metastable state. In some discussions in the literature it is not clear how this pro- 
jection is performed on rigorous grounds, although some of its physical consequences 
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are reasonable. There is indeed some criticism againt this type of approach |]TT|, ||T2f . 
It seems that more fundamental microscopic approach to deal with the tunneling in 
dissipative medium is needed. 

With this background in mind, our aim in the present work is to clarify dynamical 
aspects of the tunneling in medium: how the barrier penetration basic to quantum 
tunneling proceeds with real time. The key idea is separation of a small system from 
a larger environment, and we would like to determine the reduced density matrix for 
the small system by integrating out environment variables . This makes it easy to 
compute the penetrating flux factor for an energy eigenstate of the small subsystem. 
Our approach does not use the Euclidean technique, which in our opinion obscures 
dynamics of the time evolution. We neither assume that the tunneling system is in 
thermal equilibrium with environment, although we can discuss this case using our 
fundamental formula. Moreover, we are able to deal with both the quantum and the 
thermally activated regions in a unified way. 

The model we use to extract exact results for the barrier penetration is the 
inverted harmonic oscillator (IVHO). Since the form of the potential we use for 
exact results is very special, we cannot discuss the tunneling for general cases with 
full confidence. Nevertheless we believe that the method employed in the present 
work, especially the integral transform of the Wigner function, should be useful to 
derive approximate, yet valuable results for general potential in the semiclassical 
approach. We hope to return to a general tunneling potential in our subsequent 
work. 

We take for the environment an infinitely many harmonic oscillators of arbitrary 
spectrum. This is a standard one studied by many people in the field. The system 
of a normal harmonic oscillaor coupled to this environment is analytically solvable, 
as discussed in [fH}] , |L5[] . Our barrier penetration model corresponds to a case 



of imaginary frequency. In the present work we shall employ and extend some 
techniques we developed in the case of normal oscillator. 

The main result of the present work is summarized by a formula for the barrier 
penetration; 

f(t)\T(E)\\ /W = -^77T = ^^ ln l^)|. (1) 

This formula is applied to an eigenstate of energy E taken for the initial IVHO 
subsystem. The well known quantum mechanical penetration formula is modified by 
the environment effect in two ways; first, via the change of the original curvature uq 
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(bare parameter) to the renormalized, pole curvature u>b', 

2 1 1 

\ T ( E )\ = i + e 2-K\E\/u * I + e 2-K\E\/u B ' ( 2 ) 

with E the energy measured from the top of the potential barrier. This effect is 
essentially similar to, but numerically different from, the curvature renormalization 
effect much emphasized by Caldeira and Leggett || . In their work two cases with and 
without the friction term, but both including the curvature renormalization given 
by Ur (which is larger than ujb) are compared. The result of Caldeira and Leggett 
is understood when one writes ujb in terms of Ur. 

The second environment effect is the time dependent factor f(t) which is a ratio 
of the momentum po(t) = qo(t) to the rescaled coordinate trajectory qo(t). This 
trajectory function q (t) obyes the homogeneous Langevin equation, eq. ([43D , under 
thermal environment, being characterized by the initial energy corresponding to 
the top of the potential barrier. When the dynamical function f(t) = 1, the IVHO 
subsystem has the energy of the barrier top, and its deviation from unity is a measure 
of energy flow from the environment. When f(t) > 1, namely |po(^)| > \Qo(p)\j 
the IVHO system is excited by an energy inflow from the environment. On the 
other hand, when fit) < 1, the system is deexcited by an energy outflow. The 
main new factor f(t) deviates from unity within time range of order 1/ujb, and both 
for t <C 1/ujb and for t ^> 1/ujb the effect is small; fit) ~ 1. We find interesting 
examples in which this dynamical function exceeds unity, thus implying enhanced 
penetration, albeit for a short time of period of order 1/ujb- 

The rest of this paper is organized as follows. In Section II we explain how we 
model the environment and its interaction with a quantum system, and introduce 
the influence functional. Quantum Langevin equation is also briefly touched upon. 
In Section III we work out exact consequences for the inverted harmonic oscillator, 
and give the barrier penetration factor taking an energy eigenstate for the initial 
state. Our general result includes an integral transform of the Wigner function from 
the initial to the final one, as explained in Appendix A, . The fundamental formula 
([!]) is derived along with an explicit form of qo(t)/q (0), eq.([77|). The Ohmic or 
the local approximation (the inherently non-local term of environment interaction 
in the full Langevin equation being replaced by a few expansion terms) is shown 
to lead to some anomalous behavior of the dynamical factor. In Section IV some 
applications to the mixed initial state are discussed using the exact result of the 
inverted harmonic oscillator. Our understanding of the result of ref.f| is also made 



4 



in terms of the diagonalized curvature parameter ojb- Three appendices explain 
technical points somewhat off the main stream of arguments in the text. Appendix 
A gives the interesting form of the integral transform of the Wigner function, while 
Appendix B gives the differential form of the Fokker-Planck equation, derived both 
for the harmonic model of environment. Appendix C explains parameters necessary 
for our numerical computation of the dynamical function. 



II Model of environment and influence functional method 

We expect that the behavior of a small system immersed in thermal environment 
is insensitive to detailed modeling of the environment and the form of its interaction 
to the system. Only global quantities such as the environment temperature, the fric- 
tion, the threshold of the environment spectrum, are expected to be important. Since 
the pioneering work of Feynman- Vernon [O and Caldeira-Leggett || the standard 
model uses an infinite set of harmonic oscillators (its coordinate variable denoted by 
Q{uj)) for the environment and a bilinear interaction with the small system (denoted 
by q); 

Lq = — / dw ( Q 2 (uj) - lu 2 Q 2 (uj) J , L hlt = -q / du c{uj)Q{uj) . (3) 

We assume existence of a threshold uj c > 0. The coupling strength to the environment 
is specified by c(uj). In the present section we do not assume any special form for 
the potential of the q— system V(q). 



It is now appropriate to explain the influence functional method |I3| . The influ- 
ence functional denoted by jF[g(r) , q ,/ ( r )] results after integration of the environment 
variable Q(u) when one computes the density matrix of the entire system. Since the 
density matrix is a product of the transition amplitude and its conjugate, the path 
integral formula resulting from the environment integration has a functional depen- 
dence both on the system path q{r) and its conjugate path q'{r). We thus define, 
when an initial environment is in a mixed state given by a density matrix Pi(Qi , Q[), 
the influence functional 

T[q{T),q'{T)] = J VQ(r) J VQ'(r) J dQ, J dQ\ 

J dQ f J dQ' f 8(Q f - Q' f ) ■ K ( q(r) , Q{r) ) K* ( q'(r) , Q\r) ) p^Q, , Q' t ) , (4) 
K ( q(r) , Q(r) ) = exp ( iS Q [Q] + iS int [q , Q] ) , (5) 

S Q [Q] + Sint [q,Q)= f dT ( L Q [Q) + L im [q, , Q] ) , (6) 
Jo 



Pi(Qi,Qi') = J2 w ni>l{Qi')i> n {Qi) , (7) 

n 

where w n is the probability of finding a pure state n in the initial environment. The 
fact that there is a delta function 5{Qf — Q'f) for the environment variable at the final 
time t indicates that one does not observe the environment part at t. Throughout 
this discussion we consider a definite time interval, < r < t. 

For the thermal ensemble of a Q harmonic oscillator of frequency uj (f3 = 1/T 
being the inverse temperature), the density matrix is 



Pp{Q,Q') 



\ 1/2 

uj N 



7T coth(/3o;/2) 



exp 



UJ 



((Q 2 + Q' 2 ) cosh(/3cu) - 2QQ' 



, (8) 



2 smh((3uj) 

and one can explicitly perform the Q(uj) path integration in eq.@, since the Q(uj) 
integration is Gaussian. The result is a nonlocal action; 

n q(r) , q\r) ] = exp[ - /* dr f ds ( Z(f)a R (T - s)£(s) + ^(r)a/(r - s)X(s) ) ] , 

Jo Jo 

(9) 

with £(t) = q( T ) - q'(r) , X(r) = q(r) + g'(r) , (10) 

1 r 00 

a(r) = a R (r) + ia/(r) = — / du;a(u;) e"^ T , (11) 

27T ./— oo 

a(U )=, /; c» M ( ^_; 2+i0+ - ^ v 2 - «■> ) ■ (i2) 

Both of an and are real functions. A complex combination of these, the kernel 
function a{r) that appears in the exponent of the influence functional is the real- 
time thermal Green's function for a collection of harmonic oscillators Q{uj), which 
makes clear a relation to thermal field theory. 

The reduced density matrix p^ for the q system is defined as follows. For 
simplicity we take for the initial system a pure quantum state given by a wave 
function ip(qi); 

P iR) (<lf,<l'f) = 

Vq(r) [ Vq'(r) [ dq t [ dq[ VH) ?(t) , </(r) ] e «W«H*.M ^ft) • (13) 



Here S g [q) is the action for the q— system. This density matrix p^ may be used to 
compute physical quantities of one's interest. 

It is sometimes convenient to introduce the Wigner function by 

$>{x,p) = £* p wL + i, x -i\e-**. (14) 



( Fi) 

We shall later mention the master equation for fw ■ Here we quote for comparison a 
master equation for the Wigner function when the entire system is in a pure quantum 

state; 

dfw df w 1 / ih d ih d \ 

■P-z — ^ \ v ( x + ^r^-)- v ( x -Tr^-) ) fw- (15) 



dt dx ih \ 2 dp 2 dp 

It takes a form of infinite dimensional differential equation. 

The master equation is simplified in the semiclassical limit of h — > to 

dfw df w dV dfw nR , 



dt dx dx dp 

A great virtue of the Wigner function is that this semiclassical equation coincides 
with the Liouville equation for the distribution function of a classical statistical 
system in the phase space (x,p). It is thus easy to write down the simiclassical 
solution in the form of an integral transform, 



0) 



OO 



f w (x,p)= / dxi dp i f^(x i ,p i )6(x-x c i(t)) 5{p-p c i(t)) , (17) 



where the classical deterministic flow, (xj ,Pi) — > (x c \ ,p c i) , is defined by the classical 
mapping satisfying 

dx cl dp d fdV\ 

Pcl " ' -dT = -{-dx- _ • (18) 

\ / X=X C \ 

Although it is not our main tool of analysis, it might be of some use to recall 
the quantum Langevin equation for the model of eq.@ |l(J. By eliminating the 
environment variable Q(uj ,t) one derives the operator equation for the q variable; 

|f + ^ + 2/ o dsa I (t-s)q(s)=F Q (t), (19) 

r°° ( P(u 0) \ 

F Q (t) = - J du c{uo) ( Q(u , 0) cos(cut) + K J ; sin(c^) J , (20) 

(F Q (r)F Q ( S yy m ) env = H dur(u) cos^(r - s) coth(^) , (21) 

Jlu c Z 

where r(u) = c 2 (uj) / (2u) and ( ) env is the thermal average over the environment vari- 
ables. Thus, the kernel function a/(r) describes a nonlocal action from the environ- 
ment, while Fq(t) is a random force from the environment. The local approximation 
to ai, taking the form of a/(r) = 5uj 2 S(t) + i]S'(t) , gives 

d 2 q dV r 9 dq „ . . . . 

_*+ +fc . 8 + I ,_| = F<j(t) . (22) 
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The parameter 5uo 2 in this case is the frequency shift due to the presence of the 
environment, and 77 is the Ohmic friction. We call this approximation the Ohmic 
approximation. Perhaps more suitably, it might better be called the local approxi- 
mation. 

On the other hand, the real part of the kernel function an describes fluctuation, 
and it is related to the dissipation aj by the fluctuation-dissipation theorem; for 
their Fourier components 

- / ,x 1 „ r , ./ 2 ^(^) tanh(/V/2) 
where V denotes the principal part of integration. 



Ill Exact results for inverted harmonic oscillator 
IIIA Formalism 

We specialize the system dynamics of barrier penetration to that of the inverted 
harmonic oscillator (IVHO) given by the Lagrangian density, 

L q = ±q 2 -V q (q), V q (q) = -±u 2 q 2 , u 2 > . (24) 

There are similarities to the case of a normal oscillator of u 2 < 0, and we can take 
over some of the results derived for the unstable ( |o^o 1 2 > u 2 ) |14j], |L5] or the stable 
(0 < |co>o| 2 < 00 c) harmonic oscillator. 

The Gaussian nature of the system is a great simplification, as demonstrated 
by eq.(|]), and one may write an effective action for the q system including the 
environment effect; 

iS eS = i f dr{ l -iX + ^-iX) 
Jo 2 2 

- f dr f T ds (£(T)a R (T - s)e(s) + 1 i{r) ai {T - s)X{s)) . (25) 
Jo Jo 

The obvious linearity in the X(t) variable here gives a trivial X(r) path integration 
in the form of a delta-functional; it determines the £(r) path as 

- cu 2 £(r) + 2 jf* ds ^(s) aj(s - r) = . (26) 

The final £— path integration then leads to, as the exponent factor, 

- f dr f T dsZ(T)a R (T-s)Z(s), (27) 
Jo Jo 
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N f°° 
g{r) = — smh(ujB T ) + 2 / duj H (uj) sm(ujr) , (28) 



plus a surface term resulting after the X(t) partial integration. The function £(t) 
here is the solution of eq. (|26|) and the result of path integral must be written with a 
specified boundary condition, £(0) = £j , £(£)=£/• 

The standard technique of solving this type of integro-differential equation for 
£(t) is the Laplace transform \TA\ . We shall only summarize the final result. Two 
fundamental solutions to this equation are g(t — r) and its time derivative g(t — r) 
given by 

N 

U) B 

N=l-2 / duwH(w)<l. (29) 

JU) C 

The important properties are that g{r) is odd and g{r) is even with g(0) = 
and g(0) = 1 which gives the relation fixing N. In terms of this basic function 
g{r) a general solution to the integro-differential equation with the given boundary 
condition is 

{(TH ^ f - + ( 9(t - r) -f 9(t - r) ) {/ ' (30) 

The weight function H here is a discontinuity of some analytic function (F(z) of 
a complex variable z = uj) across the branch-point singularity along the real axis at 
u > u c , and is given by 

H(u>)= „ 2 , » 2 , rH = ^, (31) 

(uj 2 + — ll(a;)) 2 + (%r(uj)) 2 2uj 

„. . _ f 00 , , 2uj' 'r(uj') , , 

n( W ) = P / da/ , (32) 

where the integral for IT stands for its principal part. The value ujb in eq . fl28|) is 
determined as a solution for the isolated pole on the real axis at ui 2 = — ui 2 B < ; 

-u 2 B + ul + r duj ^M = 0. (33) 

U) B + U) z 



In general, one has 

<4>"o- (34) 
Since many workers in this field use a renormalized potential according to ref . M , 
we also introduce these; 

V q + V q Q = (q) + (q,Q), (35) 

Vt a) (?) = V - g 2 r du ^ = - \uj 2 R q 2 , (36) 

Jlu c UJ L 

V q { Q n \q , Q) = V qQ (q , Q) + q 2 j°° duj ^ = V qQ (q ,Q) + \ 5uj 2 q 2 . (37) 

Juj c UJ A 



We shall redefine the q— system potential using Vj Ten \q). This renormalized poten- 
tial gives an inverted harmonic oscillator with the curvature parameter lor. This 
renormalized curvature differs from the pole location of the spectral function H(u), 
namely ujb by a factor of order of the interaction coupling. In the weak coupling 
limit these two quantities are given by 



oo 



Ur = ul + 2 / du T -^- , (3£ 

UJ 



,2 2 



to r{uj) 



+ (39) 
The equation for uj 2 r is a precise definition, while the equation for uj 2 b is an approxi- 



mate relation valid for the weak coupling, the exact relation being given by eq.(33). 
In general, one can prove that 

"b<«>r, (40) 
beyond the weak coupling approximation. As an example, in the Ohmic model which 
will be discussed shortly, the relation becomes 

u B ~u R --r), (41) 



where rj is the Ohmic friction. A more precise relation in this case is eq.(|S6|). The 
relation to the bare quantity ujq is 

cu R m luq + k r] , k « . (42) 

nu 

In the infinite cutoff limit, Q —>■ oo , the quantity k is divergent. 

In Fig.l we show schematically the analytic structure of the function F(z 2 ), 
basic to the determination of the discontinuity function H(u). Unlike the case of 
the normal harmonic oscillator for which the pole location may have an imaginary 
part, the pole at z 2 = —uj 2 b appears exactly on the real axis. The other singularity 
is the branch cut starting from the threshold z 2 = uj 2 . 

The physical meaning of the basic function g{r) is better understood by solving 
the operator Langevin equation for this system; 

^-oo 2 q + 2 J t Q dra I {t-r)q{r) = F Q {r). (43) 

The quantity u 2 here should be understood as a function of u 2 B eliminating u 2 with 
eq.(^). The homogeneous solution to this Langevin equation is given by using g(t) 
and g(t); with the initial data of g(0) , g(0) = p(0) , 

q(t) = q(0)g(t)+p(0)g(t). (44) 
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The main term g(t) ~ iV smh(uBt)/uB and g(t) ~ iV cosh(ujBt) describes an 
average motion (q(t)) under the renormalized, inverted harmonic oscillator modified 
by environment, for which the original parameter \uJq\ is replaced by the new shifted 
u B . Correction to this classical motion given by the second term in eq.(28) describes 
a damped oscillation with an amplitude decreasing by an inverse power of time at 
large times . 

A straightforward calculation of X— and £— path integration finally gives an 
effective action valid for any initial state of the q— system. We first define new 
functions by 

h(u,t)= t drg{r)e-^\ (45) 
Jo 



k(u,t)= dr gM e~ luJT = tuoh(uj ,t) + g(t) e~ luJt . (46) 
J o 

With the normalization fixed by unitarity, one has for the effective action J defined 
by p( R ) = tr ( J Pi ) , 

J(X f ,t f i Xi,ti ;t) = ^je iS *, (47) 



where S c \ is given by JT 



iS* = ~^e f -^--W^^ f + i -X f i f - i -X l i, (48) 
U = {^) 2 h + I 2 -2 g -h, V=\h, (49) 

g g g 2 



W = ll a -±l 1 = l g v, (50) 
9 9 2 



h = / duo coth^-r(cj) \h(tu,t)\ 2 , (51) 

h = / duo coth^-r(cj) \k(u ,t)\ 2 , (52) 

h = / du coth^-r(uj)n[h{u,t)k*{u,t)), (53) 

{(r) = -e,^-{ / f«( t -r)-«*^W). (54) 

g(t) V 9{t) J 

For the discussion of the barrier penetration factor, we take a pure initial state 
given by a wave function if)(x) . The Wigner function in this case is 



f£>(x }P )= / -±-^^( Xi --^( Xi + -^) e - A , (55) 

J Z7TV G Z Z 

A = ~2?T + i(y9 ' Jl + y J ^( x ~ gX{ ^ + ^ J2 + 9 J 3){P- 9 Xi)) 2 

11 



+ \ ( Mx - 9Xif + J 2 (p- 9 Xif + 2J 3 (x - ^)(p- 5 Xi) ) , (56) 



/3w 



roc 

C— duo coth — r{uj)\gh{uj , t) — gk(u , t)\ 2 , 

^2.1 



hh-ir 



j, = - 



hh - n 



(57) 
(58) 



'3 -^2 — ^3 

Although it is not used in the calculation of the flux factor in the next subsection, 
it is of great theoretical interest to express our result as a transformation of the 
initial Wigner function f$ to the final one fffi- We give this in Appendix A. This 
mapping f$ — > fffi is an integrated form describing dynamics of the q— system. Its 
differential form is known as the Fokker-Planck equation, and we shall explain this 
briefly in Appendix B. 

IIIB Barrier penetration factor 

The flux at position x is computed from the formula 

dp 



I(x 



,t)= / 



to give 



I(x,t) 



dxid^i ^ 



&\ (9 



(59) 



2ixg 2 2 \g g 



■ exp 



V 



gV 



(x - gxi) - 



2g*V 



[x - gxi) 



(60) 



We use the WKB formula for energy eigenstates of IVHO. Considering the inci- 
dent left mover at x < with the unit flux, we take as the wave funtion at x > rr* 
(the right turning point) 



-^=2= exp (i f X dx'p(x')) , p(x) = J2E + uj 2 b x 2 , (61) 



where = \J2\E\/ub and T(E) is the transmission coefficient as a function of 
the energy E in a pure quantum state. This choice of the wave function gives the 
trasmission coefficient, 

TOI 2 ^ 1 + e -L/. fl - ( 62 ) 
A point which becomes important when we compare our result with those of 
other papers is how one prepares the initial state. In many past works a thermal 
equilibrium between the subsystem and the environment is often assumed, and in 
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this context it is natural to take for our choice of the pure state the reference system 
characterized by the curvature u>b, the exact pole curvature. The choice of the 
WKB wave function, using the curvature parameter u>b, fits with this picture. But 
it should be kept in mind that we may in principle take any reference curvature 
(hereafter denoted by Cj) and in these cases we replace ojb below by uj. 

Expanding wave functions in £j, we derive a general formula for the dynamical 
factor. This involves an infinite series of the expansion in £ of the initial density 
matrix element, 



n=0 



The first few terms of this series are 



Ao(xi) = \tp(xi)\ 2 , 

o L 

Computation of the flux factor is then given by 

oo 

I( X ,t) = J «(M) , 



(63) 

(64) 
(65) 

(66) 



n=0 



In(x, t)=j ^M*i)& (~ X +($- + * (- - §l 6 



2ng 



J] 



9 



r 



x exp 



9 



2 

n 



x - gxi 



(67) 



x 



/ W 



+ ia) 



da 
2^g~ 



An 



X 

j 9 



-a e 



-a 2 /2 



a + i 



h gh 
g g 2 



f3\e 



-/3 2 /2 



9 Vh, 

The second equality follows by a trivial scale change of integration variables Xi , £j. 
One may use expansion of A n {xjg + \fl\ajg) in powers of coupling, 



An I . 



-a 



An 



X 



+a: 



a g 
g 



« + ••• + 



. Aik-l) 



X 



k-1 



j/ "\gj 9 (fc-l)! " \g, 

For calculation of the penetration factor one only needs x 
function. From the WKB formula; 



9 



-a 



+ 



.(68) 



oo limit of the flux 



An\Xi 



\T(E)\ 2 d n 
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exp i 



Xi-&/2 



dx' p(x') 



,(69) 



one has 

M*) — mm 2 + o(xr 3 )} . (70) 

Only the term containing (" y j^ a ) remains here. One may then derive 

40m) — \nm-ir f^V" 1 ( 9 4 - A , (n > i) . m) 



This along with 



/(oo,o = m^)r/(o, (72) 

gives a general formula for the dynamical factor, 

/(() = X_fw_ 1 )f (_^r'= S^Stf. (73) 

This is the main result of the present work. 

The salient feature of this result is factorization; the main suppression factor 
given by \T(E)\ 2 is affected by presence of the thermal environment only via renor- 
malization effect, as will be more fully discussed shortly. The other prefactor f(t) 
carries dynamical information of the time evolution. 

The basic functions g(t) and g(t) that appear in the dynamical function f(t) are 
the solution to the homogeneous part (Fq = 0) of the Langevin equation fl43|) with 
the initial condition, g(0) = 0,g(0) = 1 . More conveniently, one may rewrite the 
dynamical function as 

/(*) = ^L Po(t)=qo(t)- (74) 
Here qo(t) is the solution with the initial condition, 

g(0) = WBg(0). (75) 
This condition corresponds to the zero energy at time t — 0, since 

H q = l -q 2 -\^ 2 B q 2 = ^- (76) 

Namely, the particle sits on the top of the renormalized potential barrier. Note 
however that the exact pole curvature ujb appears here instead of the renormalized 
ujr. The ratio of the momentum to the coordinate for f(t) is a measure of the energy 
flow from the thermal environment. Thus, the case of /(£) > (<)1 corresponds to 
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an energy inflow (outflow). Both at t <C 1/uj b and t ^> l/u B this function f(t) is 
nearly unity and it deviates appreciably from unity only within the time range of 
1/Vb- 

An explicit formula useful for detailed analysis of the dynamical function f(t) is 
q (t) = g{t)+uo B g(t) = Ne" Bt + 2 / duo H{u) ( uo cos cut + uo B sin cut ) . (77) 

Ju) c 

The first term represents a simple classical motion under the potential, modified by 
the curvature renormalization, while the second term is a further deviation due to 
the environment interaction. The environment effect for f(t) appears in two ways; 
the first via the definition of uj b determined by the potential renormalization due 
to the environment interaction. The second one is the continuous part of spectral 
integral in eq . fl77j) , and its associated deviation of the normalization factor N from 
unity. 

Before we go on to specific models of the environment, it is appropriate to discuss 
some general results. First, both in the weak coupling region and in the asymptotic 
late time the dynamical function f(t) behaves as 

2 r°° 

f(t) « /asym(t) , /asym(t) = 1 77 e~ UBt / duO H (u) (oJ 2 + 0)%) Smut . (78) 

CUijiV Ju) c 

This form has the correct asymptotic behavior, /(oo) = 1 as well as the correct 
initial behavior /(0) = 1 if the integral above is convergent. 

On the other hand, one can prove for the initial-time behavior of f(t); 

/(t)=l- WB (l-^M)t + ---, (79) 
a (0) 2 r°° 

1 V- = — duJ H(u)w ( uj 2 + u B ) > . (80) 

0J B UJ B Ju) c 

The last inequality means that for very small t , f(t) < 1, giving rise to suppression 
for barrier penetration. In order to have a convergent integral (|80D for the formula 
of 9 (0) a physical cutoff Q of the environment spectrum is necessary, thus the 
simple-minded Ohmic model without the cutoff should be treated with caution. 

The asymptotic formula (|78|) , as will be shown later, describes well, even numer- 
ically, the dynamical function at all times, including the small time. For instance, 
the expansion in time t of eq. (|78|) gives almost the same result as eq. (|79|) , except the 
missing factor N. 

We would like to stress that two properties, f(t) < 1 with /(0) = 1 for small t 
and /(oo) = 1 , are generic features for any model having a finite physical cutoff of 
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the spectral weight r(oS). As shown later, however, some approximations or models 
with the infinite cutoff violate these general properties. 

Although our fundamental formula (|72|) is derived for energy eigenstates, one 
may compute the barrier penetration factor for any mixed state by suitably weight- 
ing this rate for energy eigenstate. A salient feature of our flux formula (|7^) is that 
the dependence on the initial (system) state is given by the well known quantum 
mechanical barrier penetration factor |T(.E) 2 |, the other factor f(t) being indepen- 
dent of the initial state. This property of factorization is specific to the harmonic 
barrier. For a more general potential barrier the dynamical factor may depend on 
the energy of the initial state like f(t;E). 

We would like to stress that our result extends the result of f| in several ways. 
First, we derived the tunneling rate for any energy eigenstate, while the authors of 
ref. deals with the zero temperature limit of the mixed state in complete equilib- 
rium. Our method is also completely different from the Caldeira-Leggett's Euclidean 
approach, and our method makes it possible to discuss dynamics of the time evo- 
lution. The third point is that we derived the prefactor rigorously unlike previous 
approximate calculations. In the next section we explicitly show how we effectively 
obtain the result of Caldeira and Leggett. 

IIIC Some examples 

We would like to compute this dynanical function for a few typical examples of 
the environment spectrum. The spectral function r(u) is taken as 



(1) 


r i 




V 

= -u, 

7T 


(81) 


(2) 








(82) 


(3) 






= ^(u-u c e(u))e(\u\-u e )e(n-\u\), 

7T 


(83) 


(4) 






= ^ 2 e( W )^- M). 

7T 


(84) 



As usual, 9{x) is the step function and e(x) the signature function. The first one 
ro{uj) is the Ohmic model (without physical cutoff of the environment spectrum), 
the second ro(u;) the Drude model, and the fourth ts{uj) a super-Ohmic model, 
while the third one tt(oj) has a threshold at uj c . In Fig. 2 these spectral weights are 
schematically depicted. In the last three cases Q acts as a cutoff of the environment 
spectrum. In Appendix C we give parameters necessary to compute the dynamical 
function. 
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We note here that the Ohmic model defined here by the spectral weight ro{oS) 
should, strictly speaking, be taken as some limit of the infinite cutoff, This cutoff 
could be a straightforward frequecncy cutoff like |u;| < £1, or a smoother one as in 
the case of the Drude model of a large Q. The way how the cutoff is introduced 
does not matter provided a cutoff is there, but in some evaluation of the integral one 
should first introduce the cutoff and then take the infinite cutoff limit. 

The dynamical factor f(t) is plotted in Fig. 3—7 for these four cases. As the 
friction r\ becomes large, deviation of f(t) from unity becomes appreciable, but only 
in the time range of order 1/ub- For physical reasons we always take f] < ujb- We 
have found an interesting behavior of /(£); some models can give enhancement over 
the usual quantum formula in the time range t m 1/ub- These are the threshold 
model with a large u c and the super-Ohmic model, for which the dynamical factor 
can exceed unity. The super-Ohmic model has also a peculiar feature that the 
dynamical function f(t) can even become negative for a short time interval. The use 
of the asymptotic or the weak coupling expression for /asym(^), eq.fl75p, is compared 
to the exact result in these figures. Except at small initial times this approximate 
form gives a reasonable fit to more exact results. Presence of the normalization 
factor N in the formula is important to get a good agreement. 

The initial time behavior of the Ohmic model fl3T|) does not reflect the necessary 
condition of /(0) = 1, since f(t) at very small t is singular having no smooth deriva- 
tive at t = (the left and the right derivatives do not meet) in the case of the infinite 
cutoff limit. This implies that /(0) < 1 in this case indicates anomaly associated 
with the infinite cutoff, and should not be taken too seriously. On the other hand, 
the Drude model for a large, but finite cutoff has an expected decrease of f(t) from 
unity at small times. The larger the cutoff is, the more abrupt this decrease is, and 
a local minimum of f(t) appears at a time in proportion to 

Finally, we note how the dynamical function behaves in the Ohmic, or the local 
approximation. We wish to distinguish this Ohmic approximation from the Ohmic 
model we just discussed. First, write eq. (^2|) for IVHO; 

W +r, ti- U * q = F °®- (85) 

There is a problem of how to interpret the zero energy solution since ujb we need 
for this is not well defined. One choice might be to use the relation obtained from 
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to{oj) in the Ohmic model; 



«>B = M+ J £-j. (86) 



This relation is derived by using two exact definitions, the one for ujg (P^| ) and 
another for Ur (|38|) , along with the form of the weight r (uj). In this case the 
dynamical function is trivial; f(t) = 1. But it is by no means obvious that this 
choice is unique, since without a specific form of the weight function there is no way 
to locate the pole curvature ujr- 

Another choice is to take a more phenomenological view limiting to the local 
Langevin equation, and define the zero energy condition by taking ujr for ujr] namely, 
at time 0, 

H q = ^q 2 ~^W = 0- (87) 
The zero energy solution that initially sits on the potential top is then 

q (t) = N ( (u+ + ujr) e""* + (u- - lor) e-"+* ) . (88) 

The parameters here are given by 



w ± = ya&+^±!, (89) 
both of which are larger than ujr. It is then easy to get 

f(t) = + "r) ^ ~ ^+(^- - ur) m 

oor((uj + + ujr) + (lu_ - u R ) e~» +* ) ' 

This function behaves reasonably at initial times, having /(0) = 1. More precisely, 

f(t)^l- V t. (91) 

Its asymptotic late time behavior is given by 

/(fW(oo)(l + Ce-<"+ +"->'), (92) 
/(oo) = ^<l, c= ("++"-)("«-"-) >0 , (93) 

UJR U- [U+ + Ur) 



The function f(t) of ( p0|) is plotted in Fig.8. The property /(oo) < 1 shows an 
anomalous behavior of this local approximation, which we regard as an defect inher- 
ent in the local approximation. 



IV Applications 
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To illustrate advantages of our approach, we take up two applications of our basic 
formula, eq.([72|) along with (|73|). The first example is computation of the tunnel- 



ing rate for the kind of potential depicted in Fig. 9. Qualitatively, a new normal 
harmonic oscillator is added in the left region of the previous inverted harmonic 
oscillator. The simplest example of this class of potential is a cubic form, 

Here Vo is the barrier height measured from the left bottom of the potential, and 
ujb is taken real and positive. The important quantity added, a new frequency cu* 
in the left harmonic well, is equal to ujb in this cubic potential. Below we generally 
distinguish the two, ujb and u*, having a more general tunneling potential in mind. 

The problem we set up is to compute the tunneling rate trapped in the metastable 
state in the left oscillator part at temperature T, the same temperature as the 
environment. The setting of this problem is the same as in the Euclidean approach 



of Grabert et al.|I7[]. We assume that the temperature T < Vq with a further 
condition of Vq ^> oj*. Under this circumstance the n-th energy eigenstate of the left 
oscillator is distributed with the probability w n of the Boltzmann-Gibbs ensemble, 



e 



- n /3oj» 



We then convolute with this weight the barrier penetration factor, to get the tun- 
neling probability T, 

OO 1 

T(T , t) = (l - e-**) £ e~ n ^ \T( -V + ^(n + -))\ 2 f(t) . (96) 

We first discuss the infinite time limit, in which f(t) — > 1 . This probability has 
the zero temperature limit, 

r(o.°°)«|r(-vt + ^)|»= (97) 



which is valid for T <C u*. On the other hand, for T ^> cu* the discrete sum 

E— - . 

n=o 1 + e " 2 



Ri j 00 P -Pu* (n+l-Vo/w.) 

TiT , oo) = 2 Si „h ^ Y. '^. WMa • («» 



holds. This approximately reduces to 

\ m/r f°° dx 1 
T T , oo « 2— e"^ / - 



^w* 1 + xW/^s 

(99) 

^ BC " W ~e-<* (100) 



2 sm(Pu B /2) 
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The last formula holds for T ^> 00 b- This is the expected classical formula for 
the barrier penetration at finite temperature. More quantitatively we numerically 
computed eq.(^) to compare with various approximate formulas. The quantity 
r(T,oo), computed from eq.(|96"D with f(t) = 1, is plotted in Fig. 10. In this figure 
an approximate formula, 

^sinh(/W2) m 1 (1Q1) 

1 ' ] 00, sin(^a; B /2) 1 + e~ **(»*/fr> B )-v /uj B ) ' ^ iUi ^ 

is compared to the exact result in the case of ujb = w*- This interpolation formula 
is a simple sum of the improved high and the improved low temperature limit. 

We now present our understanding of the result of Caldeira and Leggett for the 
cubic potential V(x) (p4|). For this we extend our result of IVHO to this class of 
potential by a new formula, 

J(oo,t) =expf-2 f 2 dx J~2 (V(x) ^E) ) — — In \q (t) | , (102) 
V J xi J ujb at 

where are turning points for the energy E = — Vq solution. We take as the 
reference curvature £0 = ujb- The trajectory function qo(t) is taken as the classical, 
zero energy solution, but we ignore the dynamical function since it is almost unity 
in the equilibrium circumstance of . Thus we find for the tunneling probability to 
be compared 

I(oo,t)*\T(E)\ 2 = ew(-^ — ) ■ (103) 

V 5 uj B / 

The authors of M write the tunneling probability in terms of the renormalized 
curvature ujr and the friction 77. For the Ohmic model of small friction, we get (using 
eq.(||) with 77 <ti 0Jr), 

with 

A73^^ = 3.6^. (105) 

If we had used the IVHO potential instead of the cubic form, we would have obtained 
a numerical factor slightly different from this value, 3.6 — > n ~ 3.14. On the other 
hand, Caldeira and Leggett obtain, using a different method, 

A73 CL ^1^M^^ 6.28045^. (106) 

Thus our result gives a result numerically different by a factor of about 2 in the weak 
coupling limit. 
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On the other hand, in the large coupling limit, 77 3> 2ur, these two are 
ar 36 77V0 r/V^ 

AS = — 2~, A j B C l = 37t — . (107) 

b u R u R 

(Although our notation here suggests that this is a correction term to the main term 
of e~ 27rV °/ u>R , ^ j s actually a leading term for the case of 77 3> 2ojr.) Our result is 
about 0.8 times the Caldeira-Leggett value. 

In general, when one writes the deviation of the penetration factor from the case 
of no environment as e~ AB , one has the form, 

A73 = $(^n^. (108) 



2u R J uj r 

Both in our case and in Caldeira and Leggett the function $(a) is slowly varying 
albeit numerically different, and the deviation factor is dominated by tjVq/u r . This 
appears to be the most important dependence of parameters, the numerical factor 
being a secondary effect. 

In our interpretation of the result of Caldeira and Leggett it is crucial to use 
the pole curvature ojb as our reference, and this choice is reasonable, because it 
corresponds to an equilibrium in the zero temperature limit considered in ||. With 
the dynamical function taken as f(t) — 1, this is the only way how the friction (77) 
dependence can appear in our approach, namely via the parameter uo R in the initial 
density matrix. It is important that our inequality uj b < u R implies the general result 
of suppressed rate of tunneling in medium, the main point stressed by Caldeira and 
Leggett. We also note that detailed comparison between our result and that of || 
is possible only by assuming the relation fl86|) , specific to the Ohmic spectrum. For 
different models the difference between the two might be larger. 

In future we wish to examine further this discrepancy after we develop formalism 
for more general potentials. 

The second application is the problem of time evolution for the same potential 
as above. For this we consider the temperature range of cu* T Vo such that the 
average energy eigenstate (of E « T) may be treated semiclassically. We thus regard 
a particle in the n-th energy level (n ^> 1) of the left harmonic oscillator as moving 
almost classically with the periodic motion of frequency a;*. Each time this particle 
hits against the right barrier, it has a prescribed probability ([72] ) of tunneling into 
the far right region. Starting at time t = 0, one counts the k-th encounter with the 
barrier at time t = 2nk/uj* until the final time tf such that there are roughly tjuj^/2'K 
times of possibilities of the barrier penetration. We find it reasonable to use a reset 
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time for each encounter for the time t of eg. (|72|) . Summing up all these possibilities, 
one gets the total tunneling probability, 

r„(t / ) = |T(E n )| 2 ^7, (109) 

where / is some sort of representative value for the dyanamical factor for each en- 
counter, perhaps some average of f(t) fl75|) over one period of oscillation under the 
harmonic well like 

7 = r r t * = 7r- ( n o) 

t* JO Z7T 

Another choice for / is the dynamical function at some particular value of time 
during one period of oscillation, for instance at the classical turning point. 
For u B >• , / ~ 1, and 

r„(t)^|T(£ n )| 2 ^. (in) 

Z7T 

Thus, the total probability oc time t, and one may define the tunneling rate per unit 
time, 

V -^ = ^\T{E)\\ (112) 

This is nothing but the classic Kramers formula ||. 

On the other hand, for ujb ~ there may be a large deviation from the quantum 
mechanical formula. We plot in Fig. 11 some examples of the factor / ( 110|) as a 
function of the average time t*. In most cases studied, / is some fraction of unity, 
typically larger than 0.5. 
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Appendix A Integral transform of the Wigner function 



After some algebra, we obtain 



1 



1 



2tt j hh- n J ~™ 



dpi fw i x i J Pi) 



■ exp 



hip - Pc\) 2 + h(x - x c \) 2 - 2I 3 (x - x c i)(p - p d 



2(hh 

%ci = gxi + gpi , Pd = 9 x { + gpi . 



(113) 
(114) 



The definitions of Jj are given in the text, eqs. (|5T|) — (|53|) . The time dependent 
functions, x c \(t) , p c i{t), are homogeneous solutions to the Langevin equation ( |43| ) 
with F Q = 0. 

One may view the mapping from fjy to fffi as a kind of fluid flow. Compared 



to the classical mapping given in Section II, the quantum solution ( |113| ) in thermal 
medium is not deterministic with a broadening given by the coefficient matrix, 



{I) 




(115) 



Moreover, the initial distribution fy) itself is broadened by quantum mechanical 
effects. The peak point of the distribution is at (x c i(t) ,p c i(t)) ■ One might imagine 
that the mapping (x, ,pi) — > (x c i(t) ,p c i(t)) is not invertible due to dissipative effects 
from the environment. This is not true; the mapping is actually invertible and 



,1, i 



gxd - gpd 
g 2 -g'9 



Pi 



- 9 x d + gp c i 

g 2 -g9 



(116) 



with g 2 - g g^ 0. 



Appendix B Fokker-Planck equation 

One may derive the master equation for the reduced density matrix as described 



in ref. |L8||. Our formula for the Wigner function (|113|) is considered as an explicit 
and convenient solution to this type of master equation. From the master equation 
one can derive a Fokker-Planck equation for the Wigner function; 
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ff. AR) 

n» (t ) = -£p*4 = *lh.f? (s (ii8) 
c(o = - jfz^i (HQ) 

= if^^ ^^-^^ ^-l-^), (120) 
2 \ g9 -g 2 2g g9 -g 2 2 J 



Dxp (t) = u- 9 w + y »-™*jy w , (121) 



g 2 9 -2gg 9 +g 3 
g'g-g 2 

where coefficient functions, f2 2 (t) , C(t) , D pp (t) , D xp (t) are local functions of time 
and are written in terms of g(t) , C/ , V , W. 

Quantities that appear in this equation are well understood by writing a set of 
moment equations of low orders; 

d(x) dip) 2 



dt 
d(xp) 



dt 



= (P), ^ = -n 2 (t)(x)-C(t)(p), (122) 

= 2(xp), (123) 

= -2Q 2 (t)(xp)-2C(t)(p 2 )-4D pp (t), (124) 

= (p 2 )-n 2 (t)(x 2 )-C(t)(xp)+2D xp (t). (125) 



dt ' dt 

d{x 2 ) 

dt 
d(p 2 ) 



For instance, the quantity Q 2 (t) here is a time dependent curvature parameter mod- 
ified from the original u 2 to that in thermal medium, while C(t) is a time dependent 
friction. In similar fashons one understands D pp and D xp as fluctuations. Physi- 
cal behaviors of the harmonic oscillator system under thermal environment are all 
determined by these four quantities which are functions of the local time t. 
Limiting values relevant to large times, t ^> 1/lob, are 

tt 2 (t) ^ -lu 2 b , C(*)w0, (126) 
D„(t) « - ^ gW « ^ jH rfc cosh(^) lf(u;) , (127) 

D xp (t) & U - gW & duo cosh(— ) (u 2 + -u 2 B )H(u) . (128) 

JUlr. 2 4 



Appendix C Parameters in specific environment models 
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We give various parameters in four specific models of environment given in 
Section IIIC. These parameters are used to calculate g(t) according to 

N r°° 
g(t) = — swh(u B t) +2 dujH(uj) sin(cjt) , (129) 

H ^ = (2^2 n^L T 2 ( ,2 » ( 13 °) 
u 2 R = u 2 B + C{u B ). (131) 

The parameters are given as follows. 
(1) Ohmic model 

N= { 1 + &r) 1 ' C M=Vuj b , n(w) = 0, (132) 



g(t) = (e Wfl * - e -(-+^) . (133) 



2u 



B 



(2) Drude model 

N - + < 134 > 

^ i]uj B d t]uj 2 Q 

CM = WTn> nM = ^T^' (135) 



TV . -i / \ rfil 



2 



g(t) = — smh(u> B t) + 



ujb " v ' (a% - uJ B )(a 2 _ - uj 2 B )(a 2 + - a 2 



a 1 ) e ~"Bt + (a 2 _ ^2 ^-a+t + ^2 _ a 2_) e -«-*) , (136) 

= ^ ± /BL^I, ( 137 ) 



(3) Threshold model 



-1 



AT = l_-L_ ^ + _!_ arc tan — arctan — , (138) 

\ 7T ll 2 + UJ B TXUJb UJb KIjJb ^B J 

2tjub ( ^ ^ c 



C(u, B ) = ^ (arctan ii- arctan^) 



uj b oj B/ 
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(4) Super-Ohmic model 

= ^iin^i, nM = -^!in^!. (142) 

In calculation of g(t) for the Ohmic model one needs a frequency cutoff in in- 
termediate steps of integration, but the final result does not depend on this cutoff 
factor. For the model having a threshold and for the suer-Ohmic model we cannot 
get analytic forms of the basic function g(t). 
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Figure caption 



Fig.l 

Analytic structure giving the spectrum. The pole at — uj 2 b moves as indicated 
when the friction becomes large with fixed. 

Fig.2 

Schematic form of model spectral weights. 
Fig.3 

Dynamical function for the Ohmic model. Values of the friction relative to the 
system curvature are given for each case. Dotted lines are calculated using the 
approximate, asymptotic formula /asym(^), eq.( [78|) , in the text. 

Fig.4 

Dynamical function for the Drude model. Values of the friction and the cutoff 
relative to the system curvature are given for each case. Dotted lines are calculated 
using the approximate, asymptotic formula / asym (f), eq.(|78|), in the text. 

Fig.5 

Dynamical function for the threshold model of low threshold. Values of the 
friction, the cutoff, and the threshold relative to the system curvature are given for 
each case. Dotted lines are calculated using the approximate, asymptotic formula 
,/asym {t} j eq.(||), in the text. 

Fig.6 

The same as in Fig.5, for the threshold model of high threshold. 
Fig.7 

Dynamical function for the super-Ohmic model. Values of the friction and the 
cutoff are given for each case. Dotted lines are calculated using the approximate, 
asymptotic formula fasym(t), eq.([78p, in the text. 

Fig.8 

Dynamical function in the local, Ohmic approximation. Values of the friction are 
given for each case. 
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Fig.9 

Schematic form of a general potential. 
Fig.10 

Thermally averaged tunneling probability based on eq.fl96|) with f(t) = l in the 
text taking cu* = ojb- Crossed points are calculated using the approximate formula 

Fig. 11 

Time averaged dynamical factor for a few models. 
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